The missing link between ultraluminous X-ray sources and 

metallicity 



Michela Mapelli 

INAF-Osservatorio Astronomico di Padova, Vicolo dell'Osservatorio 5, 1-35122, Padova, Italy 

michela . mapelliOoapd . inaf . it 



ABSTRACT 

The nature of ultraluminous X-ray sources (ULXs) is still debated. Recent studies show that 
metal-poor massive stars can collapse into massive stellar black holes (MSBHs), that is black 
holes with mass > 25 M Q . Such MSBHs are sufficiently massive to explain most ULXs without 
requiring substantial violations of the Eddington limit. The recent finding of an anti-correlation 
between metallicity of the environment and number of ULXs per galaxy supports this hypothesis. 
We present the results of recent N— body simulations, including metallicity dependent stellar 
evolution, and we discuss the main pathways to produce X-ray binaries powered by MSBHs. 

Subject headings: black hole physics - stars: binaries: general - galaxies: star clusters: general - X-rays: 
binaries - methods: numerical - stars: kinematics and dynamics. 



1. Introduction 

Ultraluminous X-ray sources (ULXs) are point- 
like off-nuclear sources with X-ray luminosity, as- 
sumed isotropic, Lx > 10 39 erg s _1 . This corre- 
sponds to the Eddington limit of a > 7 M© black 
hole (BH). Various explanations have been pro- 
posed for the nature of ULXs. ULXs could be as- 
sociated with high-mass X-ray binaries (HMXBs) 
powered by stellar-mass BHs with anisotropic X- 
ray emission (e.g. King et al. 2001) or with super- 
Eddington accretion rate/luminosity (e.g. Begel- 
man 2002; Poutanen et al. 2007). Some ULXs 
(e.g. Farrell et al. 2009) could be even as- 
sociated with binaries powered by intermediate- 
mass BHs (IMBHs), i.e. BHs with mass 10 2 < 
mBH/M Q < 10 5 . Recent models (e.g., Mapelli, 
Colpi & Zampieri 2009; Zampieri & Roberts 2009) 
suggest that a fraction of ULXs may be powered 
by massive stellar BHs (MSBHs), i.e. BHs in the 
mass range 25 < ttibh/Mq < 80. These MSBHs 
might form from the collapse of massive stars, pro- 
vided that the metallicity of the progenitor star is 
sufficiently low (e.g. Belczynski et al. 2010, here- 
after B10). Finally, the objects classified as ULXs 
might actually be a mixed bag of different sources. 



The ULXs match the correlation between X- 
ray luminosity and star formation rate (SFR) re- 
ported by various studies (Grimm, Gilfanov & 
Sunyaev 2003; Ranalli, Comastri & Setti 2003). 
A number of studies suggest an (anti-)correlation 
between ULXs and metallicity, and propose that 
this may be connected with the influence of metal- 
licity on the evolution of massive stars (Pakull & 
Mirioni 2002; Zampieri et al. 2004; Soria et al. 
2005; Swartz, Soria k Tennant 2008; Mapelli et 
al. 2009; Zampieri & Roberts 2009; Mapelli et al. 
2010, hereafter M10; Mapelli et al. 2011, hereafter 
Mil). 

In this paper, we discuss the observational hints 
for a connection between ULXs and metallicity 
(Section 2), and we present the results of N— body 
simulations studying the population of X-ray bina- 
ries at various metallicities (Section 3). In Section 
4, we summarize our main conclusions. 

2. ULXs and metallicity 

Mil analyze a sample of 66 late- type nearby 
galaxies. 64 galaxies were taken from the sam- 
ple listed in Table 1 of M10. The remain- 
ing two are I Zw 18 and the interacting pair 
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SBS 0335052/SBS 0335052W1. These two objects 
are important, because they are the only galaxies 
with metallicity Z < 0.03 Z© and with published 
X-ray observations (Thuan et al. 2004). For each 
galaxy in the sample, Mil derive a fiducial value 
for the SFR, for the metallicity (adopting the 
calibration by Pilyugin & Thuan 2005; when a 
metallicity gradient is available, M10 and Mil 
adopt the value of Z at 0.7 Holmberg radii), and 
estimate the number of ULXs, Nulx, after sub- 
tracting the background contamination (see M10 
for details) . Fig. Q] shows Nulx versus SFR for 
the galaxies in Mil. This Figure points out that 
Nulx scales almost linearly with the SFR, and it 
also shows that a galaxy with relatively low metal- 
licity (Z < 0.2 Z©, black filled circles) hosts on 
average more ULXs than a galaxy with the same 
SFR but higher metallicity (Z > 0.2 Z©, red open 
circles). 

This trend is confirmed in Fig. [21 where 
Nulx /SFR versus Z is shown for the same sample 
of galaxies. The red line in Fig. [2] comes from one 
of the models described in M10. To obtain the red 
line, we assumed that Nulx = e Nbh, where Nbh 
is the number of MSBHs per galaxy, while e is the 
fraction of MSBHs that power ULXs. Nbh was es- 
timated as described in eq. 1 of M10 (adopting the 
B10 model for BH formation and a Kroupa initial 
mass function), and we assume e = 3.8 x 10 -4 , 
which matches the data of Mil. 

3. Model and simulations 

As shown in Fig. [31 the model in which MSBHs 
power (a fraction of) ULXs is able to reproduce the 
metallicity trend in observed ULXs. On the other 
hand, Linden et al. (2010, hereafter L10) perform 
population synthesis simulations of binaries with 
different metallicities, and find that the contribu- 
tion of MSBHs to bright X-ray binaries is marginal 
with respect to that of stellar-mass BHs emitting 
at super-Eddington rate. L10 explain the excess of 
ULXs in low-metallicity environments as the effect 
of metal-dependent stellar evolution, rather than 
as a consequence of the remnant mass. 

L10 analysis considers only field binaries, which 
evolve in isolation. On the other hand, most 
stars, and especially massive stars, are expected 
to form in star clusters (SCs, e.g. Lada & Lada 
2003). Furthermore, a large fraction of ULXs are 
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Fig. 1. — Nulx versus SFR in the sample by Mil. 
Red open circles: host galaxies with Z > 0.2 Z©; 
black filled circles: host galaxies with Z < 0.2 Z©. 
The solid black line is the power-law fit for the 
entire sample. The dashed red line is the power- 
law fit obtained assuming that the index of the 
power law is = 1. The error bars on both the x 
and the y axis are 1 a errors. 




Fig. 2.— Filled black circles: iV ULX /SFR versus 
metallicity Z in the sample by Mil. The error 
bars on both the x and the y axis are 1 a errors. 
The solid red line is the prediction by M10. See 
the text for more details. 
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Fig. 3. — Mass of the BH versus zero-age main 
sequence (ZAMS) mass of the progenitor star. Red 
solid line: Z = 0.01 Z©; black dotted line: Z = 
0.1 Z©; green dot-dashed line: Z = 0.3 Z©; blue 
dashed line: Z = 1 Z©. 
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associated with OB associations and with young 
intermediate-mass SCs (e.g., Goad et al. 2002; 
Zezas et al. 2002; Liu, Bregman & Seitzer 2004; 
Soria et al. 2005; Abolmasov et al. 2007; Swartz, 
Tennant & Soria 2009; Tao et al. 2011; Grise et 
al. 2012). 

SCs are often collisional environments, i.e. sys- 
tems whose two-body relaxation time is shorter 
than their lifetime. This implies that most bi- 
naries in SCs do not evolve unperturbed, but 
undergo a number of three-body encounters (i.e. 
close encounters with single stars, e.g. Sigurdsson 
& Phinney 1993). Three-body encounters can af- 
fect dramatically the orbital properties of binaries. 
Furthermore, dynamical exchanges can change the 
members of a binary. Thus, dynamics can alter 
completely the evolution of a binary with respect 
to the predictions by population synthesis codes. 

To check the combined influence of dynam- 
ics and metal-dependent stellar evolution on the 
population of X-ray sources in SCs, we ran a 
set of N— body simulations with the code Starlab 
(Portcgies Zwart et al. 2001). We included recipes 
for metal-dependent stellar evolution and for the 
formation of MSBHs (see Mapelli et al. 2013) in 
the public version of the code. Fig. [3] shows the 
mass of the BH as a function of the zero-age main 
sequence (ZAMS) mass of the progenitor star, as 
implemented in our simulations. 

We ran over 300 simulations of young (< 100 
Myr) intermediate-mass (3000 - 4000 M Q ) SCs, 
with metallicity Z = 0.01, 0.1 and 1 Z . The 
simulated SCs follow a King profile (with adimen- 
sional central potential Wo = 5 and core radius 
r c = 0.4 pc), adopt a Kroupa initial mass func- 
tion, and an initial binary fraction /pb = 0.1. The 
core collapse time for these systems is ~ 3 Myr, 
short with respect to the lifetime of massive stars. 

The simulations provide information about the 
formation and evolution of X-ray binaries, pow- 
ered by both Roche lobe overflow (RLO) and wind 
accretion. Fig. [4] shows the mass of the BH ver- 
sus the orbital period of RLO (crosses) and wind- 
accreting (circles) binaries. It is apparent that ac- 
creting binaries at low metallicity (Z = 0.01 — 0.1 
Zq) can be powered by MSBHs: 10 — 20 per cent 
of all MSBHs power wind-accreting systems, while 
5 — 10 per cent of all MSBHs power RLO systems. 
Interestingly, the vast majority (> 90 per cent) of 
accreting binaries powered by MSBHs underwent 



at least one dynamical exchange before starting 
the accretion. In most cases, the MSBH formed 
from a single star and then became member of a 
binary via exchange. 

4. Conclusions 

In this paper, we reported the results of 
N— body simulations of intermediate-mass (3000- 
4000 Mq) young SCs with three different metal- 
licities (Z=0.01, 0.1 and 1 Zq), including metal- 
dependent stellar evolution recipes and binary 
evolution. The aim is to check the importance 
of dynamics for the population of X-ray sources 
at different metallicities, and to estimate the prob- 
ability that a MSBH powers an X-ray binary. 

We showed that ^5 — 10 per cent of all the 
MSBHs formed at low metallicity (Z = 0.01 - 0.1 
Zq) power RLO systems. The vast majority of ac- 
creting binaries powered by MSBHs underwent at 
least one dynamical exchange before starting the 
accretion. Our results are in agreement with L10, 
as the number of accreting MSBHs in unperturbed 
primordial binaries is negligible. The key result of 
our simulations is that MSBHs are efficient in pow- 
ering X-ray binaries through dynamical evolution. 
This result indicates that MSBHs can power X- 
ray binaries in low-metallicity young SCs, and is 
very promising to explain the association of many 
ultraluminous X-ray sources with low-metallicity 
and actively star forming environments. 

We made use of the public software package Starlab 
(version 4.4.4) and of the SAPPORO library to run Star- 
lab on graphics processing units (CPUs). We acknowledge 
all the developers of Starlab, and especially its primary au- 
thors: Piet Hut, Steve McMillan, Jun Makino, and Simon 
Portegies Zwart. We thank the authors of SAPPORO, and 
in particular E. Gaburov, S. Harfst and S. Portegies Zwart. 
We acknowledge the CINECA Award N. HP10CXB7O8 
and HP10C894X7, 2011 for the availability of high per- 
formance computing resources and support. MM acknowl- 
edges financial support from INAF through grant PRIN- 
2011-1. 
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Fig. 4. — Mass of the BH versus orbital period in the simulated BH binaries. Filled circles: wind- accretion 
systems; red crosses: RLO systems (at the first RLO epoch). From left to right: 0.01 Z©, 0.1 Z©, 1 Z©. 
Each system can be identified by more than one point, when the period evolves significantly, as consequence 
of accretion, circularization, or dynamical interactions. 
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